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A Restricted Proof that the Weak Equivalence Principle 

* 

Implies the Einstein Equivalence Principle 

ALAN P. LIGH1MAN* and DAVID L. LEE* 

California Institute of Technology, Pasadena, California 91109 


ABSTRACT 

Schiff has conjectured that the Weak Equivalence Principle 
("WEP": free-fall trajectories independent of test-body composi- 
tion) implies the Einstein Equivalence Principle ("EEP" : all 
nongravitational laws of physics the same in every freely falling 
frame). This paper presents a proof of Schiff ’ s conjecture, 
restricted to: (i) test bodies made of electromagnetically inter- 
acting point particles, that fall from rest in a static, spheri- 
cally symmetric gravitational field; (ii) theories of gravity 
within a certain broad class — a class that includes almost all 
complete relativistic theories that we have found in the litera- 
ture, but with each theory truncated to contain only point par- 
ticles plus electromagnetic and gravitational fields. The proof 
shows that every "nonmetric" theory in the class (every theory 
that violates EEP) must violate WEP. • A formula is derived for 
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the magnitude of the violation. Comparison with the results of 
Eotvos-Dicke type experiments rules out various nonmetric theories, 
including those of Belinfante and Swihart and Naida and Capella , 

— theories that previously were believed to agree with all cur- 
rent experiments. It is shown that WEP is a powerful theoreti- 
cal and experimental tool for constraining the manner in which 
gravity couples to electromagnetism in gravitation theories. 
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I. INTRODUCTION 


In a previous paper'*' we have discussed the content and significance 
of Schiff's Conjecture. In brief, the Conjecture states that all theories 
of gravity which satisfy the Weak Equivalence Principle (WEP), i.e., pre- 
dict a unique composition- independent trajectory for any test body at a 
given point of spacetime and with a given initial velocity through that 
point, must satisfy the Einstein Equivalence Principle (EEP), i.e., must 
show that the nongravitational 1 laws of physics are the same in every freely 
falling frame. When specialized to "relativistic theories of gravity" 1 (as 
will be done throughout this paper), Schiff's Conjecture says that every 
theory satisfying WEP is necessarily a "metric theory." 1 Plausibility 
arguments (e.g., Refs. 1 and 2) have frequently been given for the Conjec- 
ture, but there have been few detailed calculations that bear upon its 
validity or invalidity. 

A particular method of attack - perhaps the only workable method for 
explicit calculations - is to attempt a proof of the Conjecture piece by 
piece: that is, first analyze test bodies with purely electromagnetic 
internal interactions and thereby attempt to show that particles and elec- 
tromagnetism must interact with gravity in the manner of metric theories 
(EEP) in order that WEP be satisfied; next analyze purely nuclear systems 
and attempt to show that nuclear fields must couple to gravity metrically, 

etc Unfortunately, for our purposes, nuclear interactions have not 

been given an adequate mathematical representation even in the absence of 
gravity; and the nonmetric theories known to us make no attempt to write 
down nuclear force laws. Hence our present program must end one way or 
another after the first stage. A proof of the first stage is still quite 
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significant, however. It will allow us to rule out various nonmetric theories 
in the literature. 

. I 

In order not to prejudice ourselves, the language and concepts used in 

the calculation will be those employed in standard classical field theory 
with gravity treated as just another ordinary field. In particular, we will 
not use such phrases as "curved spacetime" and will not make any coordinate 
transformations to real or pseudo-"freely falling frames." The concept of 
gravity as a metric phenomenon should be forced upon us by WEP. 

As spelled out in Sec. II, we shall take a nonquantum mechanical approach 
and shall use a particle rather than a fluid picture for the test body. 

Since the gravitation theories with which we attempt to tie in are largely 
classical theories, we feel that a classical approach is completely justi- 
fied and perhaps essential. There are two reasons why a particle approach 
has been taken: first, more often than not, classical field theories formu- 
late the interaction of gravity with matter in the form of point particles; 
second, a charged particle approach allows one to deal with the exact "gravi- 
tationally modified Maxwell equations" of a given theory, rather than with 
their smeared out averages. 

Our calculation is not the first of its type. For several particular 
theories, and at lower orders of approximation, the acceleration of electro- 
magnetic test bodies in a gravitational field has been previously calculated. 

•z !{. 

Nordtvedt and Belinfante and Swihart have both done calculations, to first 
order in the gravitational field potential and squared particle velocities; 

Nordtvedt for general metric theories, and Belinfante and Swihart for their 

5 

theory of gravity. In addition, Post has done a calculation, at Post- 
Newtonian order, of the acceleration of a confined quantity of electromagnetic 
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energy in a gravitational field. Had his calculation been carried to higher 
order it is conceivable he could have obtained part of our result: that 
6 = p. [cf. Eq. (21)]. 

Section II of this paper gives an outline of the assumptions, procedure, 
and techniques of our calculation, including the results; Sec. Ill presents 
the details. Section IV compares the predictions for WEP violation with the 
results of Eotvos-Dicke type experiments, and thereby rules out the non- 

Jl g. y Q 

metric theories of Belinfante and Swihart, ’ Capella, Naida, and Whitehead. 
Also discussed is the manner — both quantitative and qualitative — in which 
WEP is an experimental probe of the "gravitational-Maxwell equations," as 
contrasted to previously recognized experimental tests of those equations. 


II. GENERAL FRAMEWORK AND RESULTS 

In calculating the center of mass acceleration of an electromagnetic 
test body, we would like to set up a formalism which includes as many types 
of gravitation theories as possible, but which is not too complicated. In 
particular, our formalism should be able to deal with scalar, vector, tensor, 
scalar-tensor, etc. theories. 

We have found that all of these different types of theories can be put 
into a somewhat universal form when describing a static, spherically sym- 
metric (SSS) gravitational field - providing their dynamical law 1 for particle 
motion is derivable from a Lagrangian. (The restriction to SSS fields is 
certainly a limitation in principle, but it allows us to handle many dif- 
ferent theories at once; and, as discussed in Sec. IV, is not a limitation 
in practice. ) The quasiuniversal description of particles and electromagne- 
tism in an SSS field is as follows: 
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The motion of charged particles under the joint action of gravity and 
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the electromagnetic field A can be derived from the Lagrangian 

M* 


L = Si 


4 - 


”o k ( T - + ^ 


•d 


dt 


( 1 ) 


where we have used the bar above the L to indicate that L may be only a part 
of the total Lagrangian, and where the various symbols will be defined below 
The "gravitationally modified Maxwell equations" (GMM: Maxwell's equations 
in the presence of a gravitational field) are of the form 


V * (eE) = bnp 


( 2 ) 


v X + S/dt(eE) 


(3) 


Definitions of the quantities in Eqs. (l)-(3) and of other quantities that 
will be used in the calculation are given below: 

= spatial coordinates; they are nearly Cartesian when gravity 

is weak, 

t = a time coordinate associated with the static nature of the 

SSS field, nearly equal to proper time for slowly moving 

particles when gravity is weak, 

m , 5 rest mass of particle k, a constant, 
ok 

e = charge of particle k, a constant, 

K 

x,^(t) 5 world line of particle k, 
k 


x° = t, 

V, 2 : R.i v/ v,^ with 5., the 3-Kronecker delta, 

~k ij k k 

U(r) = a gravitational potential equal to M g /r, where M g is a con- 
stant ("active gravitational mass") characterizing the source 


h 



of the SSS field, and r is coordinate distance, 

[ (x - * (y - y s f * (* - * s )) 1/2 > fr ”> source of £leld 

point, 

V,V • = the usual differential operators of gravity free Euclidean 
space, 

g = VU - the gravitational acceleration to be expected if the 
theory in question were Newtonian theory, 

T,H,e,|a = functions of the gravitational potential U; functions that 
are arbitrary in this calculation but that have a specific 
form in each theory of gravity when the coordinate system 
has been suitably specified, 

A ^ 1 = components of an electromagnetic vector potential, a four- 
vector, 

(A ) 1 = A = spatial part of vector potential, 

5 " A 0 

2 = s k e k A* - ~k (t) ) 
p H 2 k e k &3 (^ - ?k (t) ) 

E = VA 0 - dA/dt 
B = V X A . 


(ba) 

(4b) 

(4c) 

(4d) 


Although in most theories the form of L in Eq. (l) is typical only of 
SSS fields, it turns out that all of the results we shall obtain hold even 
if U is an arbitrary, but time independent function of position. 

For an SSS field in a given theory, T, H, e, and \i will be particular 
functions of U (an hence of position). Here we a ssume that T, H, e, and ji . 
have been given and we seek the relations among them, if any, that are 
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required for compliance with WEP . It is clear from Eq. (l) that we have 
sacrificed general covariance of the particle Lagrangian in order to encom- 
pass a wide range of theories. 

Note that Eqs. (2)-(3) can be re- interpreted (different physics; same 
mathematical representation!) as the usual Maxwell equations for a permeable 
medium in which the free sources originate from charged particles labeled 
by k. Thus e and p play the role of "gravitationally induced dielectric and 
permeability parameters," respectively. We require that T, H, e, P all 
approach 1 as U vanishes so that the special relativistic limit is maintained. 

Given the SSS restriction, one may ask how general are Eqs. (l)-(3). 
Except in the most general (nonmetric) case of Jordan's theory, which is 
incomplete 1 in the sense that it involves unspecified processes of particle 
creation, all theories we know of which are complete enough to formulate 
the interaction of the electromagnetic field with gravity have GMM equations 
of the form of Eqs. (2)-(3). 12 In fact, the "e-p formulation" of the source- 
less Maxwell equations in metric theories has sometimes been used in calcu- 
lations. 13 The particle Lagrangian L [cf. Eq. (l)] also appears to be 

8 

fairly general, except for a class of theories discussed by Naida which 
includes the theory of Capella. 7 We treat the Capella-Naida theory on an 
individual basis in Sec. IV, using the methods developed in this section. 

We point out that it is sometimes necessary to perform a reformulation 
(same theory; new "mathematical representation") of a theory in order to 
put it into the form of Eqs. (l) - (3) (see, for example, the Belinfante- 
Swihart theory as analyzed in Ref. 14). Finally, we should emphasize that, 
even more important than the generality of Eqs. ( l)- (3 ), are the techniques 
and methods developed in this section, since they can also be applied on an 
individual basis to that handful of theories which is not included in 
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Eqs. (l)-(3). We now proceed with an outline of our calculations. 

Variation of Eq. (l) yields an expression for the acceleration of the 
k th particle, which, together with Eqs. (2) and (3) constitutes three 
coupled equations. We seek a perturbation solution. There are two obvious, 

small dimensionless quantities in which one could expand: the gravitational 

2 

potential U and the squared particle velocities y k . Since we prefer a 
result correct to all orders in the gravitational potential, we expand only 
in v, 2 and leave T, H, e, and \i as arbitrary functions of U. We do, however, 
expand these latter functions in a Taylor series about the instantaneous 
center of mass of the test body (defined below), i.e., 

T = T 0 + (g ' x)T 0 ' + ... , (5) 

where 

T' = dT/dU and T Q ' 5 (dT/dU) x _ Q . (6) 


We shall assume that the body is small enough so that second derivates 
of U make negligible contributions. Indeed, this is part of the definition 
of "test body" (Ref. l) and is a necessary and integral qualification in 
Schiff's Conjecture. 

We define the center of mass for the test body by the following 
sequence of equations: 


m k 


= tn 


Ok 


1 + F[U(x k )] 


+ 2 6 k E i e i ~ik' 


1 2 
+ 2 “ok-k 


1 + G[U(x k )] 


1 + K[U(x i )] + S[U(x k )] 


+ 0(m Q v ) 


(?) 


~ik ~i " 5 k ’ 


M 5 




X = M E, m. x. 
-cm k k~k 


( 8 ) 

(9) 
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Here F, G, K, S are again arbitrary functions of the potential U. (When- 
ever two indicies, e.g., i and k, occur in terms, in double or single sums, 
it is always assumed that i / k in the sum. ) Any credible result should be 
independent of the particular definition of the center of mass as long as 
it remains inside of the body, that is, the result should not depend on the 
specific forms of the functions F, G, K, and S. 

We now assume that at t = 0, the center of mass of the test body is 
momentarily at rest, at the origin of the coordinate system, 


(ijt-o = (Ut=o “ 0 • 


( 10 ) 


By differentiating Eq. (9) twice and combining with Eqs. (10), we obtain 
for the instantaneous center of mass acceleration 


where 


= M_1 ' S k\S k + ® k \X k + Z k Vk> - 


S k ^ <V dt 

m^ = dm^/dt , etc 


• • 


( 11 ) 


Return for a moment to the details of the expansion scheme. Our 
expansion is in the quantity 

= ( typical squared particle velocity) > v ^ 


(12a) 


The virial theorem guarantees that 


(typical charge of a particle) 


> 


V “ (typical mass) (typical separation of neighboring particles) ~ m^x fc | 

Thus, without serious error, we may treat both terms on the right-hand sides 

2 

of Eqs. (12a) and (12b) as 0(v ) when ordering the terms in the expansion. 


(12b) 
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Besides the dimensionless quantity v in which we do expand, and the 
dimensionless quantity U in which we do not expand, there is a third, less 
obvious dimensionless quantity: 


gs 


= |g| • (size of test body) i |g| |x^ 


(13) 


We shall expand in this quantity — independently of the v expansion — but, 

in practice, by examining powers of g rather than gs. 

Now if *X is to be body independent in general, it must be so for 
’ -cm 

each order in v 2 and each order in g, independently. Surprisingly, perhaps, 

2 

it will be sufficient to work to first order in v and to first order in g. 
The imposition of WEP at this order will force the dynamical equations (l)- 
(3) to take on metric form, thereby guaranteeing that EEP ( and hence WEP 
a fortiori ) is satisfied at all orders.' ! 

To first order in v 2 and g, after solving Eqs. (l)-(3) for a k and 
substitution into Eq. (ll), we find (details given in next section) 


5cm ~ ~ \ S ( V * 0 1) + * 


V 1 ) S i"0i v i 2 + 1 Si.k'UkJ 

+ “o ' 1 “ s i,k "Ik + M o' lcc i m Oi ( S • Jihi . 


where 


M 0 " 2 i m 0i * 


co 


= (^VxI'o'^nWo' 1 ! > 

E -|( t o 1/2h o’ 1)( 5 VV^o *\ W 


- “o' W 1 ) 


/ 1 , F ' T 1 / 2 H - — (l + G )T ' T 1//2 H _2 e 

V 1 + V F 0 T 0 0 e 0 2 K 0 J 0 0 0 € 0 

e— ■ 


-1 


0 = T ' T 

W 0 0 


^ik " e i e k'~ikl 


H 'H/ 1 + 2(1 + Fj 


-1 


0 0 


o' 


V - 2 (1 + G o) T 0 


* u - 1 ] 

H 0 j ’ 


(lb) 

(15a) 

(15b) 

(15c) 

(I5d) 

(I5e) 
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'~ik ~ e i e k^ 


?ik } l5ik 


.-3 


x M 

~ik 


(l5f ) 


Equation (lk) becomes much simplified when we use some gravitationally modi- 
fied virial relations (see Sec. III.C for details): 

< S i m 01 V iV + i (T 0 1/ V le o' 1) S l, k e i e k X ik” x ilc P Kkl'^ = 0(V 2gs) (16) 

where m, p refer to components of the appropriate vectors and <( )> denotes 
the usual time average. Using Eq. (16), Eq.' (lk) becomes 


O - - 1 * s 

• Wo H 0 ) < 5 ; ik' l ik' > 

- i V 1(I o' T 0 1/2 "o’ 2€ o" 1)( Vo" 1 • • (17) 

The first term of this acceleration is body independent (satisfies WEP); 
the second term depends on the body's self -electromagnetic energy; the third 
term depends on the electromagnetic energy, the shape of the body, and the 
orientation of the body with respect to the gravitational field gradient. 

Thus ^ X will always be body independent only if the second and third terms 

always vanish, i.e., 

(18a) 




2 € o' € o 


H o'/ H o - 2e o'/ e o - T o ,J (fA = 0 


H c/ T o " € o M, o 0 


(18b) 


(the other factors in the body dependent terms must be nonzero for correct 
Newtonian and special relativistic limits) or equivalently, 

<o7«o - ? (H o'/ H o - V>V (19a) 

^o = V<W • (19b) 
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Since we have not specified the initial location of our test body with 
respect to the external gravitating source, and Eqs. (19) should be satis- 
fied at any point we choose to deposit the body, the naught subscript can 
be removed from quantities in those equations, yielding, upon integration, 

e = CtH/T) 1 / 2 (20a) 

H = C _1 (H/T) 1/2 , (20b) 

where C is a constant. Since, "in the absence of gravity," 1 we must have 
€ = H = T = 1, C must also be unity. Therefore we finally obtain, as a 
necessary condition for our electromagnetic test body to fall with a com- 
position-independent acceleration: 

e = \x = (H/T) l/2 . (21) 

It is worth noting that, using heuristic arguments (see e.g., Ref. 15) 
about the electromagnetic energy content of atoms and the expression for 
the fine structure "constant" a in a dielectric medium 

a = (en) 1 / 2 e 2 / (efi) 

one can see why WEP should require constancy of the ratio (e/n). 

Comparison of Eqs. (21) and (l)-(3) with the discussion in Sec. III.E 
reveals that Eq. (21 ) is a necessary and sufficient condition for the 
dynamical equations (l)-(3) to take on the familiar metric form 

(22> 

F 0 ^ = bnj a . (23) 

;p 

In this metric form 

ds 2 = g Q pdx Q: dx P , (2!+a) 
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s oo T ’ 


(24b) 


,H (spherical coordinates turn out to be "isotropic 1 ), (24c) 


S ij " ‘ 6 ij 

; denotes the covariant derivate with respect to , 

F» - g'W, - A t 1 , 


^X,T 




J 01 ; S J e k 6^[x - z(s)](dx k a /ds)(- g) 1 ^ ds k . 


(24d) 

(24e) 


Note that all dependence on the arbitrary functions used in the center 

of mass definition, Eq. ( 7 ), has vanished by the time one reaches Eq. (17). 

k 2 

Higher order calculations [v or (gs) , for example] could only yield 
results consistent with Eq. (2l), since WEP at first order implies that 
gravity has a metric theory description (automatically satisfyi ng WEP) to 
all orders . 

Our theoretical results can be summarized by the following theorem: 
Theorem: Consider the class of gravitation theories that possesses a mathe- 
matical representation of the form of Eqs. (l)-(3). For that class, with 
each theory written in that representation, 


(WEP)«=->(Eq. 21)4=> 


the theory is metric with the metric 
given by Eqs. (24b)- (24c) 


III. DETAILS OF THE CALCULATION 
A. Single Particle Equations of Motion 

Variation of Eq. (l) with respect to the coordinates of particle k 
yields 

(HW -1 )a k + v k d(HW" 1 )/dt + W -1 ^(T - Hv^) = L.A. (x fc ) , (25) 


12 



where 


W = (T - Hv k 2 ) l/2 , 

L. A. (x, ) = Lorentz acceleration of particle k 
~ ~ ~ k 


(26a) 


= ( e k / m 0 k^[“ + * ~^k^) ‘ d / dt ~^] » ^ 26b ^ 

and all functions of U are evaluated on the particle's world line, e.g., 

H = h|u[ x^( t) ]| . Using Eqs. (5)-(6) and the discussion following Eqs. (13), 
we can write, to the order of our calculation. 


VH=H Q 'g etc. (27) 

We shall regard g as spatially constant [see discussion following Eq. (6)]. 
Equation (25) can then be written as 

Sk = I S (H o' v k 2 - T o' )H o ' 1 - Xk ( Xk • 

- V. (v. • a )HW“ 2 + (WH -1 )L. A. 

i ~ K "'K ~ ~ ~ 

Note that whenever functions like H, T, e, etc. occur in terms multiplied by 
g, we may evaluate them at naught, i.e., 

Hg -► H 0 § > 

because we work only to first order in g. 

2 

We further expand W in a power series in v and, since we are only 
working to 0(v 2 ), we can set W = T^ 2 in Eq. (28). This follows from the 
fact that L.A. ~ 0(v ) and from the explicit velocity dependence of other 
terms in Eq. (28). [It should be mentioned that when a term is considered 

O 

0(v ), it is not necessarily intended that the term is dimensionless, but 
2 

only that v (or the expression in Eq(l2b)) is a multiplicative factor in 
the term. The same applies to the notation 0(g).] 


- — (T ' 
2 V 0 


- v. 


V )w 


-•] 



By dotting v. into both sides of Eq. (28), solving for (g k ' Y k ), and 
substituting the result back into Eq. (28), we obtain 

* k = i g (H o’ v k 2 - T o ,)H o _1 + ^k ( Yk * sXW 1 ■ VV l) + (t1/2r 1} ~' 

+ 0(v k ) .+ 0(g 2 ) 


B. The Gravitationally Modified Maxwell Equations 

We must now solve Maxwell's equations and compute the quantity L.A. 
which occurs in Eq. (29). If Eqs. (4c) and (4d) are substituted into 
Eqs. (2) and (3) and one uses the gauge 

(ep) dcp/dt + v • A = 0 (30) 

the result is 

V^cp = epd 2 cp/dt 2 - 4irpe - e Ve * (Vtp + dA/dt) , (31a) 

V 2 A = epd 2 A/dt 2 - 4flpJ + (ep) 1 (v • A) v(ep) + M- (V X A) X • (31b) 

We can now do a perturbation solution of these equations by expanding simul- 

2 2 

taneously in powers of v and g, treating formally v g. 


-6 

II 

<* 

+ cp x + cp 2 + • 

• • ? 




(32a) 

*> 

II 

o 

+ A^ + A2 + • 

• • y 




(32b) 

cf 5 

II 

- 4rte" 1 p , 





(33a) 

II 

% 

€ p(S 2 cp 0 /St 2 ) 

-1 

" e o € o 

'[§ ’ 

(VCp Q + 

ay st)1 , 

(33b) 

II 

ep(S 2 cp 1 /St 2 ) 

-1 

' e 0 e 0 

’[g • 

(VtPi + 

dA./Bt)| etc. , 

(33c) 

ii 

o 

c 

cu 

t>* 

- 4npJ , 





(34a) 

v 2 Ai = 

eii(S 2 A 0 /Bt 2 ) 

+ (e|i) 0 “ 

X (v * 

’ 7 S + ^o"Vo'(v X A 0 ) 

X g etc. 


A. 

(29) 


(34b) 



(One should not confuse the perturbation order of A, A^, with the k com- 
ponent of the vector A^. ) 

The solution of these equations is far simpler if we remember from the 

2 

beginning that since the particle acceleration is required only to 0(v ) 
and 0(g) we need L.A. only to the same order. Remember also that 
a = o(v 2 ) + 0(g) whenever the solution of Eqs. (33)-(34) requires a particle 
acceleration as a source term (right-hand side of equations). 

We solve the equations for A first. Clearly, from the expression for 
J [ cf . Eq. (4a) ], 


*o ( 5k> - 


-1 


(35) 


Equation (35) gives the lowest order vector potential at particle k due to 
all other particles (i / k). Note that ^.(x^) is considered to be a constant 
with respect to the del -Amber tian operator acting on functions of x^. The 
above A^ can produce terms of the desired order in L.A. For example, 

e k d/dt A Q (x k ) = E i e.e k a^x^ |x fci | -1 + ... (36a) 

= S i e i®k S^^i^lSkil" + *" ( 36b ) 

2 2 

where we have substituted a. = g + 0(v ) + 0(g ). The indicated term in 

~i w 

2 

Eq. (36b) is bilinear in v and g is therefore acceptable. However it can 
be shown that no higher orders of A after A^ can contribute. For example, 
the second source term on the right-hand side of Eq. (34b) makes the contribu- 
tion, 

A x ~ 0(g) • A q ~ 0(g) 0(v 3 ) , 

L.A. ~ dA^/dt + v(v • A^) = O(gv^) + 0(g 2 v 2 ) 

From the expression for p [cf. Eq. (4b)], we can write down the lowest 


15 



order solution for the scalar electromagnetic potential: 

,-l 




( 37 ) 


The source term proportional to c^/St in Eq. (33b) doesn't contribute to 
our order of calculation. Now, define a " superpotential" X by the equation 


v 2 x = cp Q 


(38) 


Using X we can write Eq. (33b) as, to appropriate order, 

V \ = V 2 (eiia 2 x/St 2 ) - Mep.) •' v(S 2 x/St 2 ) - V 2 [e 0 _1 e 0 '(g • VX)] • (39) 

Using Eqs. (37) and (38), we obtain 

X(x k ) = \ S i e i e " 1( ~i^~ki' 

dx/dt = - \ 2 i e 1 (v i • x k .) e'^x.) I^f 1 - \ E^g • v.) |^ ki | (40b) 

a 2 x/at 2 = - •|s i e.(a. • x ki ) e' 1 (x i )|x k .|' 1 + 0(v U ) , (1+Oc) 

where we have carefully interpreted the partial time derivative on functions 

of x as acting on coordinates of particles labeled i with i 4 k. From 
~k 

Eq. (40c) it is clear that the second source term in Eq. (39) does not con- 
tribute and the remaining equation is trivially integrated to yield 


cp x = e|id 2 x/dt 2 - £q 1 € 0 ' (| * Vx) * 

Using Eqs. (40a) and (40c), Eq. (4l) becomes 

»i = - I >*0 E i e i ( 5i ' iki ) l5ul" 1 .- l £ o"%'Vi ( S 

and, using Eq. (29) for a^ 

'Pi = (J T o'^0 H 0 e oV) Z i S i^? * -ki^-wJ 


(41) 


,-l 


5ki } I ~ki ^ ~ 1 {k2) 


(43) 
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In the same manner as with the vector potential, one can show that 
5 p 2> etc. do not contribute to the Lorentz acceleration at the desired 

order. Using Eqs. (26b), (35), (37), (43), one obtains 


“ ^ e k/ m Ok^ E i 


ki 6 ' 1 (Si^Xkil' 3 - s^i^ki kkl 


+ 2(2 T o' p o H o € ° e 


-0 


o’) 


^ki “ ^ e k/ m 0k^ e ii 1 ~ki 1 


, m 


where w is as defined in Eq. (I5f). From Eqs. (29) and (6) we obtain the 
~kx 


relations 


Si = - i<V H o' > l - °' v > 


€(x ± ) = e 0 + (g - x t ) e 0 ' , 


(45a) 

(45b) 


which, when substituted into Eq. (44), yield 

"(§ ’ ~i^ki 1 ( e 0 e o' + 2 T o'^0 H 0 ) 


-1 , - 2 / 

e 0 e 0 


, , , , J 2ki e o 

L.A. (j k ) = EitejVok) T - 1 

1 l~ki 

+ 2(2 T 0 ,|a 0 H 0 " € 0 e o) E i 


Ukil 


ki 


+ 2 


-ki 


(46) 


C. Virial Conditions 


We now have enough information to derive some useful virial conditions. 
Substitution of the expression for L.A. [cf. Eq. (46)] into Eq. (29) reveals 


"okK )F -■ T o 1 / V £ o'Vi e k ( ’ ! ki )P ISkii'° + 0(g) ’ 


,-3 


(47) 


where p denotes a particular vector component. Multiplication of both sides 

£ 

of Eq. (47) with (x k ) yields 

nu,x,/ e a, P = ® 0 k d ( x k £v k P )/ dt: “ m 0k v k £v k P = T 0 e 0 E i e i e k^ X ki^ * k ^ ki ' 


Ok k k 


m 
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If we sum Eq. (U8) over the index k, use the antisymmetry of x fci , and take 
a time average, the result is Eq. (16). Summing Eq. (16) on l and p produces 
another useful virial relation: 


+ 1 ’oW 1 - 0 + 0(s) • 


(49) 


D. Center of Mass Acceleration 

We now have all of the necessary tools at our disposal for calculating 
the test body acceleration. We begin with Eq. (7). To the required order 


m k = m 0k 


F o’ ( § • + (1 + G)( Sk ' Ck J + I v k 2G o’ ( s • V 


+ 2 s i e i e k 


«o’ <* • v + s o’ ( 8 ' Ik 5 


Hikl 


( K (*l) + s (*k>) (Jik • XikllSikl" 

"k ■ "ok[ F o ,( l ' tv) + ( ?k ' !k )(1 + G) ] • 


In obtaining Eqs. (50)-(5l) we have, as before, used the fact that 
a ~ 0(g) + 0(v 2 * ). To be exact, Eqs. (29) and (be) show that 

5 k = - \ S^oV 1 * + SkilXkil' 3 + °(^ 2 ) 

Using Eqs. (50)-(52), the first two terms in the expression for 
[cf. Eq. ( 11 ) ] become 

.11 


"‘‘w ' 


— M -1 H _1 e _1 T F ' - ( 

2 0 £ 0 0 0 v 


1 + G o } T o’ H o 


(50) 

(51) 


(52) 


S . ,(1), . 
i, k~ki 


a ''Wk - 2 M_1 


V - i (1 + G o> W 1 


Vokhk • 8) 2i 


(53a) 

(53b) 


2 . 

Again using Eqs. (29) and (k6) to get the 0(gv ) contribution to a k 
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[cf. Eq. (52)], the third and last term contributing to X cm is 


M ~Vk5k = M_ 1 ij' I M o T o' (1 + VV 1 + I H o 1 

- i T 0 ' u + y] Vok v k 2 + 1 T i s i, k^ik 


H o’ (1 + F o ) 


(1 + F o )(I o" r o' 1 - M ' ls V F k • *> Z\ 


+ 2 V^kSlk ’ 


(54) 


where 


T 2 _ T 0 


1 5 T o X/ V 1(1 + F o )(£ o'V * I T o Vo -1 ) - I T o' H o' 1(1 + K o + V (55a) 

+ F o )H o' 1 (l V^o - «<fV> * £ o\‘ 


+ id + F „> ^'V 1 *, • 


(55b) 


with M Q , T) ik , u ik defined in Eqs. (15). 

Now, expand the expression for M * using Eqs. (7) and (8): 

-1 -1. \-lL 1 (1+G 0 } 2 l( 1 + K 0 + S oL 1 

M = M Q (1 + F q ) |^1 - 2 M (1 + F ) S k m Ok V k " 2 M (1 + F ) Z i, k^ikj 


+ 0(v ) + 0(g) . 


(56) 


With Eqs. (53)-(56), the expression for X^, Eq. (ll), becomes that given 
in Eq. (l4). Use of Eqs. (16) and (b9) then yields Eq. ( 17), and subse- 
quently Eq. (21). 

E. The "e - Formulation for Metric Theories 

In any static, spherically symmetric, locally Lorentz manifold with 
metric, one can introduce "spatially isotropic coordinates," for which 
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6 00 B oo 


g 0k ° * 


(r) , 


(57a) 

(57b) 


*u - - 8 u £(r) > 


5 [(x 1 - x s 1 ) 2 + (x 2 - */) + (x 3 - * s 3 ) 2 ] l/2 • ( 57c ) 


(For proof, see any standard textbook on general relativity.) For the 
problem at hand we can regard g Q0 and f as functions of U = M g /r rather 
than as functions of r. In such a coordinate system, the standard metric- 
theory Lagrangian for the motion of charged particles reduces to 


L = 2, I - m 


Ok 


i < v x k“V> 1/2 + r 


= s. 


J - - B 0k (s 00 - ^k 2)l/2 - e k*|i V k' 1 


V k 

dt ; 


(58) 


and the metric theory Maxwell equations read 

F^ - <- g)- 1 ^- g )^% - ^ , 


(59a) 


where 


J a s S k e k J* (dx k a /ds k ) B \x - x k )(- g )" 1 / 2 ds k 

= 2 k e J" S)" 1 ^ 2 6^(x - x k )(dx k a /dt) 


(59b) 


Here g = determinant of g^, and commas and semicolons denote partial and 
covariant differentiation respectively. Combining Eqs. (59) gives 


[gV^O 8) l/2 ],(3 - ll ”W 3 (; - x k )(<lx k a /dt) 


( 60 ) 


Equations (60), when written out for the diagonal, spatially isotropic metric 
of Eq. (57), have the " e - jj." form of Eqs. (2) and (3), with 


20 



etc. . 


(61) 

( 2 )- 

(62) 

one can define an "effective electromagnetic metric" by 


Sh 

II 

o 

o 

(63a) 

; ij * - s ij ; 

(63b) 


then the GMM equations will take on metric- theory form. In Eqs. (63) Y 
is an arbitrary function and reflects the well-known conformal invariance 
of Maxwell's equations. If, in addition to satisfying Eq. (62), the effec- 
tive metric determined by Eqs. (63) is correctly related to the functions 
appearing in the particle Lagrangian [cf. Eqs. (57)- (58)], then the entire 
theory of particles and electromagnetic fields can be consistently put into 
metric form. 

IV. CONCLUSIONS AND APPLICATIONS 
A. Theoretical Implications of the Results 

We have shown that, in a spherically symmetric gravitational field, a 
theory of gravity described by Eqs. (l)-(4) can be put into metric form 
(with respect to the dynamical equations for particles and electromagnetic 
fields) if and only if it satisfies the Weak Equivalence Principle. Equiva- 
lently, if such a theory is nonmetric then Eq. (21) will not be satisfied, 


and 


E i = F 0i > 


€ — U- - (f/s qq) 


1/2 


Conversely, for a theory with GMM equations of the form of Eqs. 
(3) and with 

€ = |J. 


21 



the acceleration of test bodies will have body-dependent contributions 
[cf. Eq. (17)], and WEP will be violated. The result has far reaching con- 
sequences if one accepts WEP as a valid principle: Having proved, from WEP, 
the metric nature of the GMM equations inside of an electromagnetic test 
body, one knows how to describe all gravitational-electromagnetic 
phenomena - e.g., the bending of light by the sun, electromagnetic radia- 
tion in a gravitational field, etc. 

There are two potential weaknesses of our calculation. First we have 
assumed a spherically symmetric gravitational field. Now, it is conceivable 
that a theory could be of "metric form" for spherically symmetric gravitational 
fields, but nonmetric in other cases. Such theories would have to be analyzed 
on an individual basis, to see whether their non-SSS fields violated WEP. 

However, we feel that such a theory would be difficult to formulate and, 
in fact, have seen no examples in the literature. In practical applications, 
one considers a particular nonmetric theory, solves the spherically symmetric 
problem, and finds that Eq. (21) is not satisfied, thus constituting a viola- 
tion of WEP at some order. Examples will be given below. 

A second possible weakness, discussed previously, is the limitation to 
the types of equation discussed in the beginning of Sec. II. However, except 
for the Naida-Capella nonmetric theory, discussed below, Eqs. (l)-(4) appear 
to be quite general among "complete" theories. (There are many theories 
which are not explicit as to the formulation of the GMM equations, and we 
must require that such theories be completed before given further considera- 
tion. ) 

Finally, we point out that WEP and Eq. (21) demand that the center of 
mass acceleration be body independent at each order in the external gravitational 
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potential U. As will be seen below, a given theory violating the WE P will 
do so at some order of U. To be more explicit, suppose that one expands 
the functions H, T, (i, e appearing in Eq. (17) in a power series in U, i.e., 


7 A 

H = 1 + 2?U + | &U + ... , (64a) 
T = 1 - 2QU + 2pU 2 + ... , (64b) 
e = 1 + e x U + e 2 U 2 + ... , (64c) 
ix = 1 + HjO + \i 2 U 2 + ... . (64d) 



(For the correct Newtonian limit, one must require that a = 1, but we leave 
a arbitrary here.) Each theory will yield certain values for the r's and 
T's. We have shown that nonmetric theories must have some of the r's or 
T's nonzero — the first nonzero T or T determines the order at which the 
theory violates WEP. 
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B. Experimental Verification of WEP and Applications 
of Our Calculations 

Thus far, our results have been completely within a theoretical context. 
We now investigate the experimental and practical applications. 

Experimental support for WEP comes from the type of experiment devel- 
oped by Eotvos in the late nineteenth century, and redesigned extensively 
by Dicke in the 1960's.^ The particular Eotvos-Dicke (ED) experiments of 

highest reported precision are the Princeton experiment of Roll, Krotkov, 

lfi 17 

and Dicke, ° and the Moscow experiment of Braginsky and Panov. These 

experiments measure the relative acceleration toward the sun of two dif- 
ferent substances (gold and aluminum in the Princeton experiment; platinum 
and aluminum in the Moscow experiment). The reported results are 


l^cm^Al ' ^-cm^uJ ^cro^Al ~ ^cm^uJ ^ 10 -11 


l<*cm>l 


(67a) 


l^cn^Al " ^crn^Ptl < 1Q -12 
l<*cir> I 


(67b) 


Our calculation involved a test body dropped in a static field. The 
following argument justifies direct comparison of our calculation with the 


results of the above experiments; 

(i) The 2k-hour component of the acceleration can easily be isolated 
so that the sun can really be considered as the sole external source of 
gravitation (see page 173 of Ref. 16). To make this more clear, if one uses 
the 2k-hour period variation to select out g gun from g gun + g earth , then 
Eq. (17) has body dependent terms of the form 


2b 



< '?cm' > ” Isun M 0 ^i, k^ik^ ^0 + r i^ U sun + U earth 


) + ...] 


SsunV^k^ik) Cr 0 +r i U 


sun 


+ » . . ] 


since U « 10 U , • 
sun earth 

(ii) The fact that the earth is rotating rather than at rest can only 
contribute inertial accelerations; in particular no relative accelerations 
between the two test bodies can be introduced in this manner. 

(iii) We have considered only electromagnetic test bodies; but we wish 
to apply our results to the actual atoms used in the experiments, atoms which 
have nuclear as well as electromagnetic interactions. Thus the complete 
equation for <^X ^ for realistic atoms has, in addition to the terms shown 
in Eq. (17), terms which involve nuclear energies. Is it possible that the 
nuclear and electromagnetic terms would cancel each other? The only mech- 
anism by which the terms could be combined and related is through the virial 
relations; yet an examination of Eq. (17) reveals that does not even occur 
in the electromagnetic portion of the virial relations. In particular, 
given the combined virial relations for both electromagnetic and nuclear 
interactions one could construct an infinity of different theories merely 
by changing p (and thus changing the body dependent terms in <C? cm ) > ) • Thus 
there is no credible mechanism by which nuclear and electromagnetic body 
dependent terms could conspire to cancel each other. The "electromagnetic 
violation" of WEP thus constitutes a lower limit to the total violation 
(allowing for possible nuclear violations). 

We can now ask to what order does Eq. (67) test the GMM equations of 
a theory. Equation (17) has the form 
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4m>~ i 


electromagnetic energy 


total mass 


F(H 0 , T n , e n , ii n , H o > T n » 6 0 ^ 


‘O' O’ >*0’ o 


0 


+ body independent term 


( 68 ) 


where F is a function of the indicated variables. Now, the largest contri- 
bution to the electromagnetic energy of the total atom certainly comes from 

the nuclear protons and for platinum or gold this amount to, using the 

,18 

semiempirical mass formula. 


[ (electromagnetic energy) /(total mass)] pt Qr Au 
For aluminum, the corresponding quantity is 


5 X 10 


-3 


I EM energy 

(A-V3 )ai 

EM energy 

1 total mass 

A1 (^' V3 )p t „ Au 

total mass 

L 



Pt 


2 X 10 


-3 


(69a) 


(69b) 


Noting that U Q has the magnitude 


XJ^ = potential of sun at earth ~ 10 


-8 


and using Eqs. (65) and (67), we see that current experimental accuracy 
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bears upon the T k and only for k =£ 1. The accuracy of the experiment 

must go up by a factor of IQ 7 to require that r 2 and TV , vanish ! ! Equations 

2 

(66) show that the experiment thus measures H, T, and e to 0(U ), but p 
only to 0(U). We expect that almost all theories will do well enough to 
have r Q = 0. 

Before continuing with direct applications to theories of the current 
experimental verification of WEP, let us return to Eq. (17) and analyze the 
specific way in which it constrains the GMM equations of a gravitation theory. 
The second body dependent term in Eq. (17) - the "directional Coulomb energy" 
term - involves the GMM equations only through the product ep. This particular 
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product is also equal to the square of the index of refraction, n , and is 
tested by light bending and time delay experiments (see, e.g., Ref. 20 for 
a discussion of these experiments — although in the context of metric theories). 
In fact, exploiting the "e - p" analogy for the GMM equations and taking the 
geometrical optics limit, one sees that the current experimental tests, with 
the exception of WEP, are sensitive only to the product ep — and only to first 
order in U of that quantity. On the other hand, the first body dependent 
term in Eq. (l 7) — the "nondirectional Coulomb energy" term — samples the 
GMM equations in a deeper manner, both qualitatively and quantitatively. 

Not only is e distinguished from p (magnetic and electric effects distin- 
guished) but also is e explored to second order in U (cf. the e Q ' ) for the 
current experimental verification of WEP. Thus WEP is revealed as a power- 
ful tool for probing the GMM equations — the most sensitive probe of those 
equations existing in 1973. 

On purely theoretical grounds one can require, as we have previously 
remarked, that the r's and T's vanish independently. However, in practical 
experimental applications, the second body dependent vector in Eq. (65) has 
some particular relation to the first for any given experiment. Since the 
nuclei of the atoms in the ED experiment are approximately spherical, 
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C. Applications to Specific Nonmetric Theories 


In this section we discuss WEP for three particular nonmetric theories. 
The Belinfante-Swihart and Whitehead theories have equations of the form of 
Eqs. (l)-(3). As an illustration of the formalism of Sec. IV. A and IV. B, 
the WEP violation is calculated explicitly in the case of the Belinfante- 
Swihart theory. The Naida-Capella theory, which is an apparently rare 
example of a theory not having a particle Lagrangian of the form of Eq. (l) 
in the SSS limit, is treated on an individual basis, using the techniques 
developed in Secs. II and III. 

1). 0 

1. Belinfante-Swihart Theory ’ 

An analysis of the Belinfante-Swihart theory in Ref. l4 reveals that 
its particle Lagrangian can be put into metric form with 


Sop - (l - 


. + h „ + 5 h h„ + 0(h 3 ) 

og 4 at 1 


(72) 


cyfi) 

where K is an arbitrary constant, h 5 Tj^h , and ^ is the M i- nkowski metric. 
The GMM equations are of "e - fJ." form [i.e., have the form of Eqs. (2)-(3)], 
with, in the SSS limit, 

r , 1 , 

(73a) 


e = 


M- = 


1 " 2 ( ' h 00 + h ll^ 


-1 


1 + 2^ h 00 + h ll^ 


(73b) 


In the SSS limit, h has the form 


c rP ’ 


00 0 


h ij 8 ij C l U ’ 


h 0k =0 * 


' ( 74a) 
(7to>) 
(74c) 
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( 75 ) 


where C^, are arbitrary constants, but with the implicit relation 

2K(3C 1 - C Q ) ' + C Q - 2 = 0 

in order to satisfy the Newtonian limit (g Q0 = - 1 + 2U + ...)• Defining 
T and H by comparison of Eqs. (72), (7*0 with Eqs. (24) and then evaluating 
the various T k and T fe , [ cf . Eqs. (64) and (66)], one finds 



r - — T = 
1 1 3 1 


- I C 0 (C 0 + C l> * 0 • 


(76a) 

(76b) 


In order to predict an amount of light bending and perihelion shift com- 
patible with experiment, one must require that and satisfy 


.9 £ -|(C 0 + C x - 2) S 1.1 , 

(77a) 

.8 s|(C 0 + 1) —1*3 . 

(77b) 


The combinations of C Q and occurring in Eqs. (77a) and (77b) correspond 

20 

to the 7 and 3 parameters, respectively, of the "PPN formalism" and the 
experimental limits indicated above are discussed in Ref. 20. 

Using Eqs. (7l) and (77), we find that the nonmetric theory of Belinfante 


and Swihart predicts 


4 x 10" 11 < 


• * - •• - 

<*cAu or Pt ^cnrAl 


If one requires the light bending and perihelion shift predictions of the 
Belinfante-Swihart theory to be same as in general relativity, Eq. (78) 
becomes 

i 


^cni^ A.u or Pt ^-cm^Al 
<*cm> 


6 X 10 
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Thus, the Belinfante-Swihart theory violates seriously both the Princeton 
and the Moscow versions of the ED experiment. 

9 

2. Whitehead 1 s Theory 

Synge analyzes only the motion of uncharged particles and the source- 
less GMM equations in Whitehead's theory: 

5 J* ( gQp dx a dx P ) 1//2 = 0 [Eq. (1.7) of Ref. 8] (80a) 

(g^g^F ) = 0 [Eq. (1.9) of Ref. 8] (80b) 

M’V y * 

F oe,r + F P7 ,c< + W = 0 lEq - (1 ’ 9) of Re£ ’ 81 ' (BOc) 


A straightforward generalization of these equations to include sources shows 


that the GMM equations have "e - p" form in the SSS limit, with 

€ = (- gQQf )" 1 , ( Qla ) 

p = f 2 (81b) 


[in the notation of Eqs. (57)]. Using Eqs. (17), (57), and (81) one can 
then show that 


<Xcn>, 


Au or Pt 


- <L>* 


A1 


<*cn> 


10 


-3 _d_ 
dU 


|fn(- 


8 oo f ) 


(82) 


so that, for experimentally acceptable values of g Q0 and f , this version of 

-3 

Whitehead's theory violates WEP at the order of 10 ! [Note that in Whitehead's 

theory the product ep is the same as in metric theories, so that the coeffi- 
cient of the second body dependent term in Eq. (17) vanishes identically. In 
some sense one can say that, with respect to the light bending and radar time 
delay experiments, Whitehead's theory is a metric theory.] 


30 



3. Naida-Capella Theory 


rj Q 

The nonmetric theory of Capella as completed by Naida has the follow- 
ing Lagrangian [cf. Eq. (2.1) of Ref. 7]: 

! ds |j (nQ£u a /) 1//2 + X h^u^Cq^u 0 ) 1//2 j - e J* A / x ^ ( 83 ) 


L = m. 


where q^ is the Minkowski metric and 


ds 5 (q^dx^x^) 1 / 2 , 


u Q = (dx^/ds) 


The GMM equations are of "e - |i" form [cf. Eq. (3.7) of Ref. 7] with 


e = 1 + X(h nr . + h ) , 


‘00 U 11 ' 
M- = [ 1 “ X(^QQ + 


-1 


Solutions to the SSS gravitational field equations yield 


h oo ‘ c o*' lu ' 

h ij = VVj 


(84a) 

(84b) 

(85a) 

(85b) 


where C Q and are arbitrary constants. Variation of Eq. (83) and use of 
Eqs. (85) gives the particle equation of motion [analogue of Eq. (29)] 

Sk = ij^o " C C' C 0 * 20 1 1 U 0 + C l V k " U 0 v k < 2C 0 C 1 + C 0 + 2C 1 'j 

- ^k ( ^k ‘ S)f C 0 + C 1 - 2C 1 (C 0 + c i )u oJ + H-[ l - U(C 0 + “l*]' • 


Using Eqs. (84)- (86), the GMM equations give 



( 87 ) 


L. A. (x. ) 

~ ~ v ~i 


- % + C \ 2) Z k e i e kl5lkl‘ 3 ilk 
+ i < m 0i )_1 lc i - “o (2c2 - C 0 C 1 )] I Z k lik 
- 1 (V 1 [c o + c - u o( 2c2 + c o c i>> E k Sik 

- ("Oi)' 1 C(1 - “oK'A*! • 5 k )|s lk r 3 S lk 


with C h C. + C^. 

Using the same center of mass formula as given in Eqs. ( 7 ) - ( 9 ) and the 
virial theorem 


< z i m Oi (v l ,a(v i )P [1 - V SC X + 2C 0» E t> kVk^lk^^ik^ l*tkl _3 > ' P + ° (g) (88) 


one finally obtains 

<?<=«> ’ S c o [ 1 + D o ( - 2C 1 + c o> 1 - I V 1(c o 2 + 5c i 2)u ol< z i, k1«> 

- V 1 I + 1 C 1 - 5 C 1 2 - C 0 2 - Vl "o ^i,kHtl> • (89) 

Now, with Eqs. (69)-(7l) we get 

K-cm^Pt or Au ^c n^ AjJ _ 10 _11 (1 + 3C. - 19C 2 - 5C 2 - 8C C ) ( 90 ) 

III i i U U 1 

The correct Newtonian and light bending results require, respectively: 

C Q = 1 , (91a) 

.9 < |(C 1 + 1) < 1.1 . ■ (91b) 

Equations (90) and (91) indicate then the relation 
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2 X 10' 10 s or Pt ~ ^ciAl s 2 x u f 10 . (92) 

<0 

Thus the Naida-Capella nonmetric theory seriously violates both the 
Princeton and Moscow versions of the ED experiment. 
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